Molecular Dynamics Simulation Study of the Selective Inhibition of Coagulation Factor IXa over Factor Xa.

Thromboembolic disorders, arising from abnormal coagulation, pose a significant risk to human life in the modern world. The FDA has recently approved several anticoagulant drugs targeting factor Xa (FXa) to manage these disorders. However, these drugs have potential side effects, leading to bleeding complications in patients. To mitigate these risks, coagulation factor IXa (FIXa) has emerged as a promising target due to its selective regulation of the intrinsic pathway. Due to the high structural and functional similarities of these coagulation factors and their inhibitor binding modes, designing a selective inhibitor specifically targeting FIXa remains a challenging task. The dynamic behavior of protein-ligand interactions and their impact on selectivity were analyzed using molecular dynamics simulation, considering the availability of potent and selective compounds for both coagulation factors and the co-crystal structures of protein-ligand complexes. Throughout the simulations, we examined ligand movements in the binding site, as well as the contact frequencies and interaction fingerprints, to gain insights into selectivity. Interaction fingerprint (IFP) analysis clearly highlights the crucial role of strong H-bond formation between the ligand and D189 and A190 in the S1 subsite for FIXa selectivity, consistent with our previous study. This dynamic analysis also reveals additional FIXa-specific interactions. Additionally, the absence of polar interactions contributes to the selectivity for FXa, as observed from the dynamic profile of interactions. A contact frequency analysis of the protein-ligand complexes provides further confirmation of the selectivity criteria for FIXa and FXa, as well as criteria for binding and activity. Moreover, a ligand movement analysis reveals key interaction dynamics that highlight the tighter binding of selective ligands to the proteins compared to non-selective and inactive ligands.

The Interaction of Factor Xa and IXa with Non-Activated Antithrombin in Michaelis Complex: Insights from Enhanced-Sampling Molecular Dynamics Simulations.

The interaction between coagulation factors Xa and IXa and the activated state of their inhibitor, antithrombin (AT),have been investigated using X-ray diffraction studies. However, only mutagenesis data are available for non-activated AT. Our aim was to propose a model based on docking and advanced-sampling molecular dynamics simulations that can reveal the conformational behavior of the systems when AT is not binding a pentasaccharide. We built the initial structure for non-activated AT-FXa and AT-FIXa complexes using HADDOCK 2.4. The conformational behavior was studied using Gaussian accelerated molecular dynamics simulations. In addition to the docked complexes, two systems based on the X-ray structures were also simulated, with and without the ligand. The simulations revealed large variability in conformation for both factors. In the docking-based complex of AT-FIXa, conformations with stable Arg150-AT interactions can exist for longer time periods but the system also has a higher tendency for reaching states with very limited interaction with the "exosite" of AT. By comparing simulations with or without the pentasaccharide, we were able to gain insights into the effects of conformational activation on the Michaelis complexes. RMSF analysis and correlation calculations for the alpha-carbon atoms revealed important details of the allosteric mechanisms. Our simulations provide atomistic models for better understanding the conformational activation mechanism of AT against its target factors.

A molecular dynamics simulation study to understand the effect of cholesterol and tissue factor palmitoylation on tissue factor-factor VIIa-factor Xa ternary complex in different lipid environments.

BACKGROUND: Tissue factor (TF), a transmembrane glycoprotein, plays a profound role in the formation of the tissue factor-factor VIIa (TF-FVIIa) complex that initiates factor Xa (FXa) generation followed by thrombin activation and clot formation. Previous reports suggest that TF-FVIIa coagulant activity at the cell surface may be affected by various processes, including changes in cholesterol content and posttranslational modifications of TF. Numerous studies were conducted but yielded inconclusive results about the effect of cholesterol on TF expression. OBJECTIVE: The present study aimed to understand how cholesterol affects structural modulations on the tissue factor-factor VIIa-factor Xa ternary complex (TF-FVIIa-FXa). Additionally, we aimed to illustrate the effect of palmitoylation on the Cys245 residue of TF and understand its structural implications on the TF-FVIIa-FXa. METHODS: We set up the following 4 systems in different lipid environments: TF-FVIIa-FXa in POPC:POPS (CS), TF-FVIIa-FXa in POPC:POPS:CHOL (CSL), Palmitoylated TF-FVIIa-FXa in POPC:POPS:CHOL (CSLP), and Palmitoylated TF-FVIIa-FXa in POPC:CHOL (CLP), respectively, and subjected them to molecular dynamics simulation. RESULTS: Hydrogen-bond and contact probability analysis were performed between various important domains of TF-FVIIa-FXa and notable novel interactions: Asn93FVIIa:L-Lys48TF, Arg178FVIIa:H-Asp95FXa:B, Lys20FVIIa:H-Glu193FXa:A, Arg178FVIIa:H-Asp97FXa:B, and Arg153FVIIa:H-Gln135FXa:B have been reported. The protein stability study implies that the CS and CLP systems are thermodynamically less stable than CSL and CSLP systems. CONCLUSION: Analysis of molecular dynamic simulation data suggests that the presence of cholesterol and palmitoylation may contribute to structural rigidity, stability, and compactness of key domains of TF-FVIIa-FXa by augmenting protein-protein and protein-lipid interactions.

Investigation of the tight binding mechanism of a new anticoagulant DD217 to factor Xa by means of molecular docking and molecular dynamics.

A new promising drug candidate DD217 has been proposed recently as a potent anticoagulant acting on factor Xa (fXa) target. It exhibits the lowest concentration of doubling the prothrombin time among the known anticoagulants. In order to explain the efficacy of DD217 in terms of molecular interactions with its target we studied the hypothesis of the tight binding mechanism by means of molecular dynamics simulations and statistical analysis of the trajectory. The conducted analysis confirms the significant contributions to the MM/GBSA estimated binding free energy of the S4 pocket residues as well the crucial role of establishing the hydrogen bonds between the ligand and the backbone amides of Gly216 and Gly218 of the target. The simulation results support the hypothesis of the tight binding mechanism of DD217 to fXa.Communicated by Ramaswamy H. Sarma.

Variability in the composition of porcine mucosal heparan sulfates.

Commercial porcine intestinal mucosal heparan sulfate (HS) is a valuable material for research into its biological functions. As it is usually produced as a side-stream of pharmaceutical heparin manufacture, its chemical composition may vary from batch to batch. We analysed the composition and structure of nine batches of HS from the same manufacturer. Statistical analysis of the disaccharide compositions placed these batches in three categories: group A had high GlcNAc and GlcNS, and low GlcN typical of HS; group B had high GlcN and GlcNS, and low GlcNAc; group C had high di- and trisulfated, and low unsulfated and monosulfated disaccharide repeats. These batches could be placed in the same categories based on their 1H NMR spectra and molecular weights. Anticoagulant and growth factor binding activities of these HS batches did not fit within these same groups but were related to the proportions of more highly sulfated disaccharide repeats.

Combinatorial assembly, traceless generation and in situ evaluation of inhibitors for therapeutically relevant serine proteases.

A combinatorial method was devised and applied for the design and identification of substrate-analogue inhibitors of therapeutically relevant serine proteases, such as thrombin and factor Xa. We conceptualized imino acid derived diketomorpholines as generally applicable key intermediates prepared through solid-phase synthesis and prone to be cleaved with primary amines in a traceless fashion. The approach led to a compound library whose members were prepared under bioassay-compatible conditions and directly subjected to the in situ evaluation, allowing a fast prediction of hit compounds. Highly active inhibitors for serine proteases of the coagulation cascade have been identified. The most potent dual inhibitor, 16K, has a binding affinity of 23.9 nM to thrombin and 32.8 nM to factor Xa.

Effect of variations in the conserved residues E371 and S359 on the structural dynamics of protein Z dependent protease inhibitor (ZPI): a molecular dynamic simulation study.

Protein Z (PZ) dependent protease inhibitor (ZPI) is a natural anticoagulant inhibiting blood coagulation proteases fXa and fXIa. Despite being a member of the serpin superfamily, it possesses unique structural features such as activation by PZ, regulating its inhibitory function. In order to understand the Reactive Centre Loop (RCL) dynamics of ZPI, which is absolutely critical for its activity, we performed Molecular Dynamics (MD) simulation on ZPI and its E371 and S359 variants located at important conserved functional sites. Unexpectedly, the RCL of E371 variants, (E371K, E371R, and E371Q), were shown to be very stable due to compensatory interactions at the proximal end of RCL. Interestingly, RCL flexibility was shown to be enhanced in the double mutant K318E-E371K due to the repulsive effect of increased negative charge on top of the breach region. Principal component analysis (PCA) coupled with residue wise interaction network analysis(RIN) revealed correlated motion between the RCL and the PZ binding regions in the WT. However, a loss of regulation in correlated motion between RCL and PZ binding hotspot Tyr240 in the double mutant was also observed. Additionally, the S359F and S359I mutations resulted in increased RCL flexibility owing to the disruption of stabilizing hydrogen bonding interaction at the distal end of strand S5A. Thus, the current study proposes that the overall stabilizing interactions of S5A is a major regulator of proper loop movement of ZPI for its activity. The results would be beneficial to engineer activity compromised ZPI as a prophylactic agent for the treatment of hemophilia.Communicated by Ramaswamy H. Sarma.

Induced forms of α2-macroglobulin neutralize heparin and direct oral anticoagulant effects.

Alpha2-macroglobulin (α2M) is a physiological macromolecule that facilitates the clearance of many proteinases, cytokines and growth factors in human. Here, we explored the effect of induced forms of α2M on anticoagulant drugs. Gla-domainless factor Xa (GDFXa) and methylamine (MA)-induced α2M were prepared and characterized by electrophoresis, immunonephelometry, chromogenic, clot waveform and rotational thromboelastometry assays. Samples from healthy volunteers and anticoagulated patients were included. In vivo neutralization of anticoagulants was evaluated in C57Bl/6JRj mouse bleeding-model. Anticoagulant binding sites on induced α2M were depicted by computer-aided energy minimization modeling. GDFXa-induced α2M neutralized dabigatran and heparins in plasma and whole blood. In mice, a single IV dose of GDFXa-induced α2M following anticoagulant administration significantly reduced blood loss and bleeding time. Being far easier to prepare, we investigated the efficacy of MA-induced α2M. It neutralized rivaroxaban, apixaban, dabigatran and heparins in spiked samples in a concentration-dependent manner and in samples from treated patients. Molecular docking analysis evidenced the ability of MA-induced α2M to bind non-covalently these compounds via some deeply buried binding sites. Induced forms of α2M have the potential to neutralize direct oral anticoagulants and heparins, and might be developed as a universal antidote in case of major bleeding or urgent surgery.

Enoxaparin Dose Requirements to Achieve Therapeutic Low-molecular-weight Heparin Anti-factor Xa Levels in Infants and Young Children.

INTRODUCTION: Enoxaparin is commonly used to treat pediatric thrombosis. Several small retrospective studies have suggested that infants and young children require higher enoxaparin doses to achieve therapeutic anti-factor Xa levels compared with adults. MATERIALS AND METHODS: This is a retrospective study of hospitalized children who received enoxaparin for the treatment of thrombosis at a free-standing children's hospital. The primary objective was to ascertain the enoxaparin dose required to achieve an anti-factor Xa level of 0.5 to 1.0 U/mL among 4 age groups in a large cohort of infants and young children between 60 days and 5 years of age. RESULTS: A total of 176 infants and children were evaluated. The majority of patients were less than 1 year of age (n=104). An inverse relationship between enoxaparin dose needed to achieve therapeutic anti-factor Xa levels and patient age was noted, particularly in the first year of life. Patients who were 60 days to less than 7 months at the time of enoxaparin initiation (n=73) required the highest mean dose among the age groups at 1.73 mg/kg subcutaneously every 12 hours (P<0.0001). CONCLUSION: Infants and young children require higher doses of enoxaparin to achieve therapeutic anti-factor Xa levels compared with adults.

The 3-O-sulfation of heparan sulfate modulates protein binding and lyase degradation.

Humans express seven heparan sulfate (HS) 3-O-sulfotransferases that differ in substrate specificity and tissue expression. Although genetic studies have indicated that 3-O-sulfated HS modulates many biological processes, ligand requirements for proteins engaging with HS modified by 3-O-sulfate (3-OS) have been difficult to determine. In particular, the context in which the 3-OS group needs to be presented for binding is largely unknown. We describe herein a modular synthetic approach that can provide structurally diverse HS oligosaccharides with and without 3-OS. The methodology was employed to prepare 27 hexasaccharides that were printed as a glycan microarray to examine ligand requirements of a wide range of HS-binding proteins. The binding selectivity of antithrombin-III (AT-III) compared well with anti-Factor Xa activity supporting robustness of the array technology. Many of the other examined HS-binding proteins required an IdoA2S-GlcNS3S6S sequon for binding but exhibited variable dependence for the 2-OS and 6-OS moieties, and a GlcA or IdoA2S residue neighboring the central GlcNS3S. The HS oligosaccharides were also examined as inhibitors of cell entry by herpes simplex virus type 1, which, surprisingly, showed a lack of dependence of 3-OS, indicating that, instead of glycoprotein D (gD), they competitively bind to gB and gC. The compounds were also used to examine substrate specificities of heparin lyases, which are enzymes used for depolymerization of HS/heparin for sequence determination and production of therapeutic heparins. It was found that cleavage by lyase II is influenced by 3-OS, while digestion by lyase I is only affected by 2-OS. Lyase III exhibited sensitivity to both 3-OS and 2-OS.

Regulation of factor V and factor V-short by TFPIα: Relationship between B-domain proteolysis and binding.

Coagulation factor V (FV) plays an anticoagulant role but serves as a procoagulant cofactor in the prothrombinase complex once activated to FVa. At the heart of these opposing effects is the proteolytic removal of its central B-domain, including conserved functional landmarks (basic region, BR; 963-1008 and acidic region 2, AR2; 1493-1537) that enforce the inactive FV procofactor state. Tissue factor pathway inhibitor α (TFPIα) has been associated with FV as well as FV-short, a physiologically relevant isoform with a shortened B-domain missing the BR. However, it is unclear which forms of FV are physiologic ligands for TFPIα. Here, we characterize the binding and regulation of FV and FV-short by TFPIα via its positively charged C-terminus (TFPIα-BR) and examine how bond cleavage in the B-domain influences these interactions. We show that FV-short is constitutively active and functions in prothrombinase like FVa. Unlike FVa, FV-short binds with high affinity (Kd ∼1 nM) to TFPIα-BR, which blocks procoagulant function unless FV-short is cleaved at Arg1545, removing AR2. Importantly, we do not observe FV binding (µM detection limit) to TFPIα. However, cleavage at Arg709 and Arg1018 displaces the FV BR, exposing AR2 and allowing TFPIα to bind via its BR. We conclude that for full-length FV, the detachment of FV BR from AR2 is necessary and sufficient for TFPIα binding and regulation. Our findings pinpoint key forms of FV, including FV-short, that act as physiologic ligands for TFPIα and establish a mechanistic framework for assessing the functional connection between these proteins.

Crystal Structure of Isoform CBd of the Basic Phospholipase A2 Subunit of Crotoxin: Description of the Structural Framework of CB for Interaction with Protein Targets.

Crotoxin, from the venom of the South American rattlesnake Crotalus durissus terrificus, is a potent heterodimeric presynaptic ß-neurotoxin that exists in individual snake venom as a mixture of isoforms of a basic phospholipase A2 (PLA2) subunit (CBa2, CBb, CBc, and CBd) and acidic subunit (CA1-4). Specific natural mutations in CB isoforms are implicated in functional differences between crotoxin isoforms. The three-dimensional structure of two individual CB isoforms (CBa2, CBc), and one isoform in a crotoxin (CA2CBb) complex, have been previously reported. This study concerns CBd, which by interaction with various protein targets exhibits many physiological or pharmacological functions. It binds with high affinity to presynaptic receptors showing neurotoxicity, but also interacts with human coagulation factor Xa (hFXa), exhibiting anticoagulant effect, and acts as a positive allosteric modulator and corrector of mutated chloride channel, cystic fibrosis transmembrane conductance regulator (CFTR), implicated in cystic fibrosis. Thus, CBd represents a novel family of agents that have potential in identifying new drug leads related to anticoagulant and anti-cystic fibrosis function. We determined here the X-ray structure of CBd and compare it with the three other natural isoforms of CB. The structural role of specific amino acid variations between CB isoforms are analyzed and the structural framework of CB for interaction with protein targets is described.

In Silico Evaluation of the Effectivity of Approved Protease Inhibitors against the Main Protease of the Novel SARS-CoV-2 Virus.

The coronavirus disease, COVID-19, caused by the novel coronavirus SARS-CoV-2, which first emerged in Wuhan, China and was made known to the World in December 2019 turned into a pandemic causing more than 126,124 deaths worldwide up to April 16th, 2020. It has 79.5% sequence identity with SARS-CoV-1 and the same strategy for host cell invasion through the ACE-2 surface protein. Since the development of novel drugs is a long-lasting process, researchers look for effective substances among drugs already approved or developed for other purposes. The 3D structure of the SARS-CoV-2 main protease was compared with the 3D structures of seven proteases, which are drug targets, and docking analysis to the SARS-CoV-2 protease structure of thirty four approved and on-trial protease inhibitors was performed. Increased 3D structural similarity between the SARS-CoV-2 main protease, the HCV protease and α-thrombin was found. According to docking analysis the most promising results were found for HCV protease, DPP-4, α-thrombin and coagulation Factor Xa known inhibitors, with several of them exhibiting estimated free binding energy lower than -8.00 kcal/mol and better prediction results than reference compounds. Since some of the compounds are well-tolerated drugs, the promising in silico results may warrant further evaluation for viral anticipation. DPP-4 inhibitors with anti-viral action may be more useful for infected patients with diabetes, while anti-coagulant treatment is proposed in severe SARS-CoV-2 induced pneumonia.

Immobilisation of an anti-platelet adhesion and anti-thrombotic drug (EP224283) on polydopamine coated vascular stent promoting anti-thrombogenic properties.

Current vascular drug-eluting stents based on immuno-proliferative drugs would reduce the rate of in-stent restenosis (ISR) but may be associated with a higher risk of acute stent thrombosis due to non-selective activity. In this paper, we aimed to develop a polydopamine (PDA) coated chromium­cobalt (CoCr) stent functionalised with EP224283 (Endotis Pharma SA), which combines both a GPIIbIIIa antagonist (tirofiban moiety) and a factor Xa inhibitor (idraparinux moiety) to reduce acute stent thrombosis. PDA-coated chromium­cobalt (CoCr) samples were first immersed in a polyethylenimine (PEI, pH 8.5) solution to increase amine function density (36.0 ± 0.1 nmol/cm2) on the CoCr surface. In a second step, avidin was grafted onto CoCr-PDA-PEI through the biotin linkage (strategy 1) or directly by coupling reactions (strategy 2). The HABA titration proved the fixation of biotin onto CoCr-PDA-PEI surface with a density of 0.74 nmol/cm2. The fixation of avidin was demonstrated by water contact angle (WCA) and surface plasmon resonance (SPR). SEM micrograph shows the flexibility of the thin layer coated onto the stent after balloon inflation. Independently of the strategy, a qualitative SEM analysis showed a reduction in platelet activation when the molecule EP224283 was immobilised on avidin. In parallel, the measurement of anticoagulant activity (anti-Xa) revealed a higher anti-factor Xa activity (2.24 IU/mL vs. 0.09 IU/mL in control) when EP224283 was immobilised on avidin. Interestingly, after seven days of degradation, the anticoagulant activity was persistent in both strategies and looked more important with the strategy 2 than in strategy 1. Throughout this work, we developed an innovative vascular stent through the immobilisation of EP224283 onto CoCr-PDA-PEI-(avidin) system, which provides a promising solution to reduce ISR and thrombosis after stent implantation.

Synthesis, Docking, and In Vitro Anticoagulant Activity Assay of Hybrid Derivatives of Pyrrolo[3,2,1-ij]Quinolin-2(1H)-one as New Inhibitors of Factor Xa and Factor XIa.

Coagulation factor Xa and factor XIa are proven to be convenient and crucial protein targets for treatment for thrombotic disorders and thereby their inhibitors can serve as effective anticoagulant drugs. In the present work, we focused on the structure-activity relationships of derivatives of pyrrolo[3,2,1-ij]quinolin-2(1H)-one and an evaluation of their activity against factor Xa and factor XIa. For this, docking-guided synthesis of nine compounds based on pyrrolo[3,2,1-ij]quinolin-2(1H)-one was carried out. For the synthesis of new hybrid hydropyrrolo[3,2,1-ij]quinolin-2(1H)-one derivatives, we used convenient structural modification of both the tetrahydro- and dihydroquinoline moiety by varying the substituents at the C6,8,9 positions. In vitro testing revealed that four derivatives were able to inhibit both coagulation factors and three compounds were selective factor XIa inhibitors. An IC50 value of 3.68 µM for was found for the best factor Xa inhibitor and 2 µM for the best factor XIa inhibitor.

Dirofilaria immitis possesses molecules with anticoagulant properties in its excretory/secretory antigens.

Dirofilaria immitis is a parasitic nematode that survives in the circulatory system of suitable hosts for many years, causing the most severe thromboembolisms when simultaneous death of adult worms occurs. The two main mechanisms responsible for thrombus formation in mammals are the activation and aggregation of platelets and the generation of fibrin through the coagulation cascade. The aim of this work was to study the anticoagulant potential of excretory/secretory antigens from D. immitis adult worms (DiES) on the coagulation cascade of the host. Anticoagulant and inhibition assays respectively showed that DiES partially alter the coagulation cascade of the host and reduce the activity of the coagulation factor Xa, a key enzyme in the coagulation process. In addition, a D. immitis protein was identified by its similarity to the homologous serpin 6 from Brugia malayi as a possible candidate to form an inhibitory complex with FXa by sodium dodecyl sulfate polyacrylamide gel electrophoresis and mass spectrometry. These results indicate that D. immitis could use the anticoagulant properties of its excretory/secretory antigens to control the formation of blood clots in its immediate intravascular habitat as a survival mechanism.

Innovative Three-Step Microwave-Promoted Synthesis of N-Propargyltetrahydroquinoline and 1,2,3-Triazole Derivatives as a Potential Factor Xa (FXa) Inhibitors: Drug Design, Synthesis, and Biological Evaluation.

The coagulation cascade is the process of the conversion of soluble fibrinogen to insoluble fibrin that terminates in production of a clot. Factor Xa (FXa) is a serine protease involved in the blood coagulation cascade. Moreover, FXa plays a vital role in the enzymatic sequence which ends with the thrombus production. Thrombosis is a common causal pathology for three widespread cardiovascular syndromes: acute coronary syndrome (ACS), venous thromboembolism (VTE), and strokes. In this research a series of N-propargyltetrahydroquinoline and 1,2,3-triazole derivatives as a potential factor Xa (FXa) inhibitor were designed, synthesized, and evaluated for their FXa inhibitor activity, cytotoxicity activity and coagulation parameters. Rational design for the desired novel molecules was performed through protein-ligand complexes selection and ligand clustering. The microwave-assisted synthetic strategy of selected compounds was carried out by using Ullmann-Goldberg, N-propargylation, Mannich addition, Friedel-Crafts, and 1,3-dipolar cycloaddition type reactions under microwave irradiation. The microwave methodology proved to be an efficient way to obtain all novel compounds in high yields (73-93%). Furthermore, a thermochemical analysis, optimization and reactivity indexes such as electronic chemical potential (µ), chemical hardness (η), and electrophilicity (ω) were performed to understand the relationship between the structure and the energetic behavior of all the series. Then, in vitro analysis showed that compounds 27, 29-31, and 34 exhibited inhibitory activity against FXa and the corresponding half maximal inhibitory concentration (IC50) values were calculated. Next, a cell viability assay in HEK293 and HepG2 cell lines, and coagulation parameters (anti FXa, Prothrombin time (PT), activated Partial Thromboplastin Time (aPTT)) of the most active novel molecules were performed to determine the corresponding cytotoxicity and possible action on clotting pathways. The obtained results suggest that compounds 27 and 29 inhibited FXa targeting through coagulation factors in the intrinsic and extrinsic pathways. However, compound 34 may target coagulation FXa mainly by the extrinsic and common pathway. Interestingly, the most active compounds in relation to the inhibition activity against FXa and coagulation parameters did not show toxicity at the performed coagulation assay concentrations. Finally, docking studies confirmed the preferential binding mode of N-propargyltetrahydroquinoline and 1,2,3-triazole derivatives inside the active site of FXa.

Elucidating the tight-binding mechanism of two oral anticoagulants to factor Xa by using induced-fit docking and molecular dynamics simulation.

Communicated by Ramaswamy H. Sarma.

Practical High-Quality Electrostatic Potential Surfaces for Drug Discovery Using a Graph-Convolutional Deep Neural Network.

Inspecting protein and ligand electrostatic potential (ESP) surfaces in order to optimize electrostatic complementarity is a key activity in drug design. These ESP surfaces need to reflect the true electrostatic nature of the molecules, which typically means time-consuming high-level quantum mechanics (QM) calculations are required. For interactive design much faster alternative methods are required. Here, we present a graph convolutional deep neural network (DNN) model, trained on ESP surfaces derived from high quality QM calculations, that generates ESP surfaces for ligands in a fraction of a second. Additionally, we describe a method for constructing fast QM-trained ESP surfaces for proteins. We show that the DNN model generates ESP surfaces that are in good agreement with QM and that the ESP values correlate well with experimental properties relevant to medicinal chemistry. We believe that these high-quality, interactive ESP surfaces form a powerful tool for driving drug discovery programs forward. The trained model and associated code are available from https://github.com/AstexUK/ESP_DNN.